Creation of a Photonic Time-bin Qubit via Parametric Interaction of Photons in a 

Driven Resonant Medium 
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A novel method of preparing a single photon in temporally-delocalized entangled modes is pro- 
posed and analyzed. We show that two single-photon pulses propagating in a driven nonabsorbing 
medium with different group velocities are temporally split under parametric interaction into well- 
separated pulses. As a consequence, the single-photon "time- bin-entangled" states are generated 
with a programmable entanglement, which is easily controlled by driving field intensity. The ex- 
perimental study of nonclassical features and nonlocality in generated states by means of balanced 
homodyne tomography is discussed. 

PACS numbers: 42.50.Dv, 03.67.-a, 03.65.Ta 



I. INTRODUCTION 



Entanglement and nonlocal correlations, besides their 
fundamental importance in the modern interpretation of 
quantum phenomena [3, 0], are the basic concepts for 
realization of quantum information procedures [3[. The 
entanglement between matter and light states is an es- 
sential element of quantum repeaters [J] , the intermedi- 
ate memory nodes in quantum communication network 
aimed at preventing the photon attenuation over long 
distances. The two-photon entanglement is a crucial in- 
gredient for quantum cryptography [a, la] , quantum tele- 
portation 0, H, 0] , and entanglement swapping [n], EH , 
which have been successfully realized during the last 
decade by utilizing two approaches, one based on contin- 
uous quadrature variables and the other using the polar- 
ization variables of quantized electromagnetic field [121 ]. 
An essential step has been recently made in this direc- 
tion by implementing robust sources producing the pairs 
of photons which are entangled in well-separated tempo- 
ral modes (time-bins) [HI]. It has been shown [l3|, [3] 
that this type of entanglement, in contrast to other ones, 
can be transferred over significantly large distances with- 
out appreciable losses, thus being much preferable for 
long-distance applications. From the fundamental view- 
point, of special interest is a single photon delocalized 
into two distinct spatial fT{| or temporal modes. In the 
last decade, the concept of single particle entanglement 
has been an object of intensive debates [H, Il7l . [l8l [l9j . 
Although, there was some criticism in literature [181 ] con- 
cerning whether a single degree of freedom can be entan- 
gled with itself, it is now well recognized that a single 
photon state delocalized spatially or temporally in the 
two modes is entangled and nonlocal 

[MtHGJl- More- 
over, one-particle entangled qubit has been used to de- 
velop and further study of quantum cryptography [2fj| . 
quantum computing with linear optics [l5| or teleporta- 



tion [2l|, [2J, |23(. The robust criteria [2J, |25| are now 
available for verification of entanglement and nonlocal- 
ity of a correlated two-mode quantum state of light via 
testing the Bell's inequality that has been recently real- 
ized experimentally [2a , l27f by performing the homodyne 
detection of delocalized single-photon Fock states and re- 
constructing the corresponding Wigner function from ho- 
modyne data (see also [28j). 

Two approaches have been hitherto developed for 
preparation of a single-photon in two distinct temporal 
modes. In first one a time-bin qubit is created with use of 
linear optics by passing a short pulse from a spontaneous 
parametric down conversion (SPDC) source through 
Mach-Zehnder interferometer with different-length arms 
[l3l ]. The second approach is based on conditional mea- 
surement on quantum system of entangled signal-idler 
pairs generated via SPDC of two consecutive pump pulses 
in a nonlinear crystal, when a detection of one idler pho- 
ton tightly projects the signal field into a single-photon 
state coherently delocalized over two temporal modes 
[27l ]. However, the both methods are confronted with 
severe challenges. The main limitation is that the light 
emitted via SPDC has too broad linewidth (~ lOnm) 
and low spectral brightness to be able to excite atomic 
species. Additionally, due to short coherence time (~ 
femtosecond) the photon waveforms are not or hardly 
resolvable by existing photodetectors, as well as their co- 
herence length is small for long distance quantum com- 
munication. 

In this paper a novel method free from the above draw- 
backs is discussed for dynamical preparation of photonic 
time-bin qubit 
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where | 0)t and | l)t denote Fock states with zero and one 
photon, respectively, at the time t and | n | 2 + | r 2 | 2 = 1. 
The basic idea is to create a parametric interaction be- 
tween two single-photon pulses, which propagate in a res- 
onantly driven medium without absorption and at low, 
but different, group velocities. Then, due to the cyclic 
parametric conversion of the fields and the group de- 
lay, each pulse experiences a temporal splitting into well- 
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separated subpulses. Moreover, since the process is com- 
pletely coherent, at the output of the medium the time- 
delocalized and entangled single-photon state is formed. 
One obvious limitation of this resonant process is that it 
is effective in a relatively narrow frequency range asso- 
ciated with specific atoms. We note, however, that re- 
cently a source of narrow-bandwidth, frequency tunable 
single photons with properties allowing exciting the nar- 
row atomic resonances has been created [HI, Ho|. Thus, 
our mechanism is a robust source for temporally entan- 
gled narrow-bandwidth single-photons. Another impor- 
tant advantage is a generation in a simple manner of 
any desired entanglement by controlling the driving field 
intensity. In section II we describe a three-level model 
parametric interaction between two quantum fields prop- 
agating in a driven medium under the conditions of elec- 
tromagnetically induced transparency (EIT). Then, in re- 
alistic approximations, we obtain an analytical solution 
for the field operators and calculate the output intensi- 
ties of the fields showing the splitting of an initial single- 
photon pulse into two well separated temporal modes. 
In section III we analyze the entanglement characterizing 
our single-photon state by employing the Bell's inequality 
proposed by Banaszek and Wodkiewicz [24[ and show an 
unambiguous correspondence of this inequality violation 
to the degree of two-mode single-photon entanglement. 
Finally, in section IV we summarize our conclusions. 



II. THREE-LEVEL MODEL OF PHOTON 
PARAMETRIC INTERACTION 

We consider an ensemble of cold atoms with level con- 
figuration depicted in Fig.l. Two quantum fields 
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£1,2(2, t) exp[i(fci i2 2 - u»i,2*)] + h.c. 



co-propagate along the z axis and interact with the atoms 
on the transitions — > 1 and — > 2, respectively, while 
the electric-dipole forbidden transition 1 — > 2 is driven 
by a classical and constant radio-frequency (rf ) field with 
real Rabi frequency £1 inducing a magnetic dipole or an 
electric quadrupole transition between the two upper lev- 
els, and V is the quantization volume taken to be equal 





FIG. 1: (a) Level scheme of atoms interacting with quantum 
fields Ei : 2 and classical rf driving field of Rabi frequency f2. 
(b) Geometry of fields propagation. 



to interaction volume. The electric fields are expressed in 
terms of the operators £i(z,t) obeying the commutation 
relations (see Appendix A) 



L 



[£i(z,t),£f(z,t')} = -5 l0 8{t-t') 



(2) 



where L is the length of the medium. We describe the 
latter using atomic operators a a p(z, t) — ^ | a )i((3 \ 

2 i=l 

averaged over the volume containing many atoms N z = 
jrdz 3> 1 around position z, where N is the total number 
of atoms. In the rotating wave picture the interaction 
Hamiltonian is given by 



H = -h^- I dz[gi£ia w e 



ik\Z 
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(3) 



Here fc|| = k^e z is the projection of the wave- vector of 

the driving field on the z axis, g a = p a oy // uJi/(2heoV) 
is the atom-field coupling constants with fi a p being the 
dipole matrix element of the atomic transition a — > /?. 
For simplicity, we discuss the case of exactly resonant 
interaction with all fields and, therefore, put in Eq.(3) 
the frequency detunings equal to zero, neglecting so the 
Dopplcr broadening, which in a cold atomic sample is 
smaller than all relaxation rates. Then, using the slowly 
varying envelope approximation, the propagation equa- 
tions for the quantum field operators take the form: 



dz + -c8i) £ ^ t)= ^-^ e 



Ft 



(4) 



(5) 



where Fi(z,t) are the commutator preserving Langevin 
operators, whose explicit form is given below. 

In the weak-field (single-photon) limit, the equa- 
tions for atomic coherences p oi = <r ai e~ lkiZ , i = 
1,2 and p\ 2 = ai 2 e~ l ( k2 ~ kl ) z are treated per- 
turbatively in £±. 2 . In the first order only (aoo) — 1 is 
different from zero and for these equations we obtain: 

d 

^7 Ah = -r 1 p Q1 +ig 1 £ 1 a 0Q + inp 02 e tAkz -ig 2 £ 2 p 21 (6) 
at 



~Qj.Po 2 = -r 2 /5 2 +ig 2 £ 2 (Too +«£lpoie lAkz -ig\£\f>\ 2 (7) 
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P12 = -I12P12 - igi£'iPo2 + «.92<?2/5io (8) 
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Here Ak = k 2 ~ ki — k» is the wave- vector mismatch and 
r 12 and r 12 are the transverse relaxation rates involving, 
apart from natural decay rates 71,2 of the excited states 1 
and 2, the dephasing rates in corresponding transitions. 
The latter are caused by atomic collisions and escape of 
atoms from the laser beam. However, in the ensemble of 
cold atoms the both effects are negligibly small compared 
to 71,2, so that 1^2=71,2/2 and r 12 =(7i + 72)/2. 

Further, we assume that the phase-matching condition 
Afc = is fulfilled in the medium. Then, the solution to 
Eqs.(6-8) to the first order in E\,i is readily found to be 

r . .r> 2 -r 2 a s n s , 2m 9- 

poi = ^fc + .-^-^-fc - -g 2 £ 2 + -^92^2, 

(9) 



/5o2=/5oi(1^2), D = n 2 +T 2 



(10) 



where, for simplicity, the optical decay rates are taken to 
be the same: Ti = T2 = T. The first terms in right hand 
side (RHS) of Eqs.(9,10) are responsible for linear absorp- 
tion of quantum fields and define the field absorption co- 
efficients Ki — g 2 TN/c£l 2 upon substituting these expres- 
sions into Eqs.(4,5). Here the condition of electromagnet- 
ically induced transparency (EIT, refs. [UUll) ^ ^ ^1,2 
is assumed to be satisfied for both transitions coupled to 
the weak-fields. Note that the three-level configurations 
0-2-1 and 0-1-2 form the A- and ladder EIT-systems, re- 
spectively, with the same decoherence time T _1 [33||. The 
second terms in RHS of Eqs.(9,10) represent the disper- 
sion contribution to the group velocities of the pulses, 
while the two rest terms describe the parametric interac- 
tion between the fields. We require that the photon ab- 
sorption be strongly reduced by imposing the condition 
KiL -C 1. Another limitation follows from AlueitT > 1 
indicating that the initial spectrum of quantum fields is 
contained within the EIT window Aujeit = ^ 2 /(^V^) 
[3^ . where T is a duration of weak-field pulses, a = NoL 
is the optical depth, a = ^ A 2 is the resonant absorption 
cross-section, and Af is the atomic number density. Fi- 
nally, the length of the pulses has to fit the length of 
the medium: Tvi < L with i?j = cfl 2 /g 2 N being the 
group velocity of the i-th field. Taking into account that 
KiL ~ r 2 a/fi 2 , this set of limitations yields 



Q 2 

r 2 



>> a 



and 



1 Tvi 
-== < — 1 < 1 



(11) 



It is worth noting that upon satisfying the conditions 
(11), the dominant contribution to the parametric cou- 
pling between the photons is the third term in RHS of 
Eq.(9,10), because in this case the last term becomes 
strongly suppressed by the factor il 2 T/T >> 1. 

It is useful at this point to consider numerical estima- 
tions. The sample is chosen to be 87 Rb vapor with the 
ground state 5Si/ 2 (F g — 2) and exited states 5P 3 / 2 (F ( 
2), 5P 3 / 2 (F e — 3) being the atomic states and 1,2 in 
Fig.l, respectively. Using the following parameters - light 
wavelength A ~ 0.8/Ltm, r = 2ir x 3 MHz, atomic density 



TV - 10 12 cur 3 in a trap of length L ~ 100 fjcm, ft ~ 10r, 
and the input pulse duration T ~ 2-^3ns, we find a ~ 16, 
v 2 ~ 10 4 m/s, vi ~ 0.3v 2 , and KiL ^ 0.1. All of the pa- 
rameters we use in our calculations appear to be within 
experimental reach, including the initial single-photon 
wave packets with a duration of several nanoseconds sat- 
isfying the narrow-line limitation discussed above. 

The noise operators Fi t % in Eqs.(4,5) have the proper- 
ties [H 

(Ft(z,t)) = (F l (z,t)F l (z,t')) = (F i (z,t)F j (z,t'))=0 



(F i (z,t)Ff(z,t')) = 2^d ij S(t-t') 
c 

showing that in the absence of photon losses the noise 
operators Fi give no contribution. Then, in this limit the 
simple propagation equations for the field operators are 
obtained: 



d_ 

dz 



j_d_ 

v\ dt 



£i(z,t) = —i/3£2 



o -1 a 

dz v 2 at' 



(12) 



(13) 



where j3 = gig2N/c$l is the parametric coupling con- 
stant. In Appendix A we show that these equations pre- 
serve the commutation relations (2). Note that for the 
parameters above, the parametric interaction between 
the photons is sufficiently strong: (3L ~ 3. 

The formal solution of Eqs.(12,13) for the field opera- 
tors in the region ^ z ^ L is written as 



£i(z,t) = £i(0, t-z/vi)+ / dx{£i(0, t-z/v 



Av jt 



dz 



i/3 £j(0, t — z/vi 



Avi 



■x) JoW}, (14) 



where i,j = 1,2 and j ^ i. The Bessel function Jo WO 
depends on z via ip = 2(3y / x/z — x) , Avij = Vi — Vj is 
the difference of group velocities. 

We are interested in the evolution of the input state 
I ipin) —\ li)<8 I O2) consisting of a single-photon wave 
packet at wi frequency, while lo 2 field is in the vacuum 
state. The similar results are clearly obtained in the case 
of one input photon at lu 2 frequency. We assume that 
initially the u\ pulse is localized around z — with a 
given temporal profile /i(t) : 

(0 I &((),*) I = (0 I &((),*) I li> = fi{t) (15) 

= where fi(t) is normalized as j- J \ fi(t) \ 2 dt = 1. In free 
space, £i(z, t) — £i(0, t — z/c) and we have 

(0\£i(0,t-z/c) I l x )=h{t-z/c). (16) 
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The intensities of the fields at any distance in the region 
$J z ^ L are given by 



(Ii(z,t)) = | (0 | £i{z,t) | Vm) 



(17) 



Using Eqs.(14-16) and recalling that (0 | £ 2 (0, i) | Vin) = 
0, we calculate (7j) numerically and show in Fig. 2 the 
output pulses at z = L for the three values of f2 and for 
Gaussian input (at z = 0) pulse /i(t) = Cexp[— 2i 2 /T 2 ], 
where C is some normalization constant. For one-photon 
initial state, as is the case here, one can clearly see that 
the second field is not practically generated, thus demon- 
strating that our scheme enables to prepare a single- 
photon in a pure temporally-delocalized state with an 
efficiency ~ 100%. Moreover, depending on the driving 
field intensity, a different degree of initial pulse splitting 
is attainable. It is easy to check that the total number 
of photons which is determined by the areas of the corre- 
sponding peaks is conserved upon propagation through 
the medium. To show this we introduce the dimensionless 
operators for numbers of photons that pass each point on 
z axis over the whole of time 



hi{z) = - 



dtEf(z,t)Ei(z,t) 



(18) 



With taking into account that (0 | £i(z,t — > ±oo) | 
V'm) = 0, the conservation low for mean photon num- 
bers results from Eqs.(12) and (13) 



d_ 

dz 



(m(z) + n 2 (z)) = 



(19) 



where i%i(z) = (ipi n \ hi{z) \ Vm)- Since in our 
case n 2 (0) = and n 2 (L) is negligibly small, we have 
n\{L) ~ ni(0) = 1. Thus, in the considered scheme 
we are able to preserve the output wi-field in a single- 
photon state while modulating its amplitude to get a de- 
sirable spatio-temporal distribution. Further, only two 
well-separated output temporal modes at w\ frequency 
are produced. To understand the physics of this splitting, 
let us discuss the structure of solution (14) for £i(z, t) in 
detail. The first term in this equation represents the lo\- 
pulse in the absence of ui 2 generation. In such a case, the 
group velocity of the pulse is slowed down to v\ << c 
under the conditions of EIT realized via ladder system 0- 
1-2. However, the input w\ photon can also be converted 
to one photon of the lu 2 field which is emitted on the 
dipolc-allowcd transition 2^0 and propagates in the 
medium at a group velocity v 2 << c established under 
EIT in the A- scheme 0-2-1. In its turn, the u> 2 photon 
is transformed back into u)\ photon, which precedes the 
signal u)\ photon owing to v 2 > V\. This process is de- 
scribed by the second term in Eq.(14). The last term in 
this equation corresponds to generation of an oj\ photon 
by an input lo 2 photon and gives no contribution in our 
case. An important point here is that an atom, being 
excited to the upper state 2 upon absorbing the initial 
u>i photon and one photon of the drive field, can return 
back to the ground state in two ways: 1) by emitting one 



photon of drive field and an u>i photon and 2) by emitting 
an lo 2 photon on the transition 2^0. The competition 
between the two processes evidently leads to a destruc- 
tive interference between the modes propagating in the 
u>i channel with different group velocities and ultimately 
gives rise to the pulse temporal splitting. In Appendix B 
we show that, indeed, the contributions of the two pro- 
cesses into lo\ photon wavefunction are of opposite signs. 
The separation of u)\ temporal modes depends clearly on 
the relative velocity of quantum fields, the larger the ra- 
tio v 2 /vi, the larger the group delay and the larger the 
separation of the two output lo\ pulses. On the contrary, 
in the limit of equal group velocities v 2 = v\ = v the 
propagating pulses experience no splitting, as it follows 
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FIG. 2: Temporal shapes of quantum fields (t is given in units 
of input pulse duration T) at the output of the medium z = L 
for three values of Q. In these figures, solid curves represent 
the uii pulses, dashed curves show the ui2 field generated in the 
medium, and dotted lines display the initial Gaussian pulse 
at loi frequency with T = 2ns propagating in the medium in 
the absence of parametric interaction (3 = 0. For the rest of 
parameters see the text. 
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from Eqs.(14), which in this case are reduced to 

£i(z, t) = £i(0, t)cos{I3z) - iij(0, r)sin(Pz) (20) 

where r = t — z/v, j =/= i. 

We finish this consideration with a short remark about 
the dynamics of frequency conversion from lj\ to u>2 and 
back to u>i in dependence on the travelled distance. For 
equal group velocities v\ = V2, a simplest result fol- 
lows from Eq.(20) with (0 | 8 2 (0, t) \ ip in ) = 0. It is 
seen that the complete conversion of the input aj± pho- 
ton to one photon of the L02 held occurs at the distance 
Zmax{v\ = V2) = 7r/(2/3), which is proportional to fl and 
thus exhibits the square-root dependence on the drive 
field intensity. An essentially different picture is ob- 
served for vi v 2 - In this case the efficient frequency 
conversion from loi to w 2 takes place up to a distance, 
where the time delay between the signal lj\ and para- 
metrically generated u>2 pulses becomes comparable to 
the pulse duration T. The numerical calculations show 
that, for the parameters above, this happens roughly at 
Zmax{v\ 7^ V2) — (v2 — V\)T ~ 20/im. At this point, the 
L02 field reaches its maximal value with conversion effi- 
ciency ~ 0.13. For z > z max (vi 7^ V2) the two pulses are 
no longer overlapped in time, and the L02 pulse is trans- 
formed to the fast lo\ pulse almost completely. Only a 
tiny part on the leading edge of the w 2 pulse leaves the 
medium with the group velocity u 2 - With further propa- 
gation in the medium, the energy is merely pumped from 
the slow uj\ pulse to the fast one, while the w 2 field re- 
mains negligibly small. This occurs as long as the fast lo\ 
pulse becomes sufficiently strong in order to generate a 
new L02 pulse. According to our analysis, the correspond- 
ing distance is approximately 130/im showing that in a 
sufficiently large interval of propagation lengths only the 
u>i field is present in the medium in the form of two well 
separated temporal modes. This result provides a wide 
choice of the length of atomic sample that is important 
for further applications of the proposed mechanism. 

It must be noted that the considered system is capable 
of fully entangling two single-photon pulses at different 
frequencies lo\ and 0J2 in the case of input state | ipi„ ) = 
I 1 1 ) <S> I I2 ) - The study of this problem is, however, 
beyond the scope of the present paper and its results will 
be published elsewhere. Here we note only that in this 
case two time-bin qubits at lo\ and 102 are generated, be- 
ing at the same time strongly correlated with each other. 
This correlation is clearly seen from the particular result 
of Eq.(20). 

III. NONLOCALITY IN GENERATED 
STATE: ANALYSIS OF BELL'S INEQUALITY 

Now we discuss whether nonlocal correlations arise be- 
tween the two generated temporal modes of a single uj\- 
photon and how they can be verified experimentally We 
first note that the entanglement is not apparent in the 
single-photon wavefunction <$>i(z, t), which is represented 



as a sum (see Eq.(B3)), rather than mixture product 
of two orthogonal mode functions <&f (z,t) and $f(z,t). 
This does not mean, however, that the entanglement is 
absent in the photon state. It can appear in second quan- 
tization formalism, where the solution Eqs.(14) for the 
field operators has been found. The point is what rep- 
resentation of Fock space has to be chosen to make this 
effect visible, since, as it has been shown in [36j, the 
entanglement with vacuum and nonlocality in a single- 
photon state is not a property of the Fock space in gen- 
eral, but appears if a specific irreducible representation is 
chosen, although the same physics follows from all Fock 
representations in the sense that the experimental test of 
entanglement and nonlocality is performed by mapping 
the quantum field state into the state of the matter (sin- 
gle trapped atoms or ions, atomic ensembles, quantum 
dots, etc.) and the results of such measurements are the 
same independent of the representation. Below we de- 
scribe the output wi-photon state in terms of quantized 
temporal modes and choose the relevant representation 
for the field quantization. 

Using the definitions of photon number operators 
Eq.(18) and wavefunction Eq.(Bl), we obtain the ui\ 
single-photon output state as 

I U)z=l = I / dt^{L,t)£+{L,t) I 0) (21) 

Normalization requires that 

{ J dt\^(L,t) | 2 =1 

This condition, the left hand side of which coincides ev- 
idently with the mean photon number ni(L), manifests 
along with ni(0) = j- J dt \ f\(t) | 2 = 1 the photon num- 
ber conservation law. Remind that for intermediate val- 
ues of z (0 < z < L) 712(2) is no longer zero and hence 
n\(z) < 1. At these distances the incoming uj\ field is 
converted, completely or partially, into LO2 optical mode 
that enables frequency conversion and redistribution of 
quantum information between different quantum fields. 
This mechanism will be discussed elsewhere. 
Let us rewrite the Eq.(21) in the form 

I li>*=£ = | / dt(<P((L,t) + <S>f(L,t))£+(L,t) I 0) 

and introduce the operators of creation of single-photon 
wave packets associated with orthogonal set of mode 
functions &\(L,t) (3?j], where i labels the members of 
the denumerably infinite set. For i — F, S these opera- 
tors are given by 

$,S = N f,I J dt^- S (L,t)£+(L,t) (23) 
with the normalization constants 

N l = ^{J dt\^\{L,t)\ 2 )- 1 (24) 
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These operators create the single-photon states in the 
usual way by operation on the vacuum state | 0) 

ct I 0) =| (25) 

and have the standard boson commutation relations 

[c h cf]=5 ij (26) 

Note that this definition of quantum temporal modes is 
only useful, if one can perform the local measurements on 
these modes such that they are spacelike separated. That 
is why the requirement for the modes be well separated 
is important. Now, for the algebra (26) we choose the 
representation of infinite product of all vacua 

io}=nio>*=ioMo)5 n ( 2? ) 

However, since in our problem we deal with two modes, 
while the other modes are not occupied by the photons 
and, hence, are not taken into account during the mea- 
surements, the vacuum may be reduced to | 0) =| 0}f \ 
0)s- Then the single-photon state (22) can be written as 

I 1i) 2 =l = nc+ I 0) F I 0) 5 + r 2 | 0) F c+ | 0) s 

= n | 1) F I 0) s + r 2 | 0) F I l)s (28) 

which is just the state (1) with r\, 2 = 

J £ J dt\ $f" s '(i,t) | 2 . The remarkable property 
of chosen representation (27) is that in this case the 
entanglement in the single-photon state is entirely 
converted into nonlocal entanglement between the 
atoms [36] and, hence, the entanglement in the field 
state reproduces adequately the expected results of any 
measurement one may perform on atomic systems. The 
amount of entanglement in the state (28) is simply 

E=-[r\ log 2 r\ + (1 - r\) log 2 (l - r 2 )] (29) 

calculated as E = — Tr$[p log 2 p] with p = Trp | |- 

It is easy to check that E is maximal E max = 1 for ri 
i 

Now we pass to discussion of nonlocal correlations in 
the state (28). Note that the single-photon states are 
completely described by their Wigner function, whose re- 
markable property is that it takes negative values around 
the origin of phase space for the complex field ampli- 
tude. The negativity of the Wigner function is the ulti- 
mate signature of non-classical nature of these states. Be- 
sides, the nonlocality of quantum correlations in single- 
photon qubit is directly evident from the violation of 
Bell's inequality formulated for two-mode Wigner func- 
tion. Specifically, we will employ the criterion proposed 
by Banaszek and Wodkiewicz [24[ allowing a Bell test 
with high levels of violation. The ability of this approach 



has recently been demonstrated in the case of SPDC tem- 
porally entangled single photon [27| • As has been shown 
in Ref-Hil, the local theories impose the bound 

- 2 ^ B ^ 2 (30) 

where the combination B has the form: 

7T 2 

B = — [W(0, 0) + W(ai,0) + W(0, a 2 ) - W(ai,a 2 )] 

Here W(ai, a 2 ) is the Wigner function of two temporal 
modes calculated for complex amplitudes a, = Xi + iyi 
with Xi and yi, i — 1,2, being the quadratures of the i-th 
mode. The Wigner function for the state (28) is obtained 
to be: 

W( ai ,a 2 ) = -5. [2 I r iai +r 2 a 2 | 2 _i] e -2|«i| 2 -2|a 2 | 2 

(31) 

which is always negative W^O, 0) < independent of r\ 
and r 2 . The amplitudes r± and r 2 depend on the driving 
field intensity and arc calculated numerically. 

Then the strongest violation of inequality (30) is 
achieved when a = a± = a 2l for which case the com- 
bination B takes the form: 

B = -l + (4J-2)e- 2J -[4J(r 1 + ^/l - r 2 ) 2 -l] e - 4J (32) 

where J — a 2 . Its behavior is plotted in Fig. 3 as a 
function of J for the values of r\ corresponding to three 
output states of lo\ photon depicted in Fig. 2. It is ev- 
ident that the maximal violation, which is about ten 
percent (2.2 compared to a classical maximum of 2), is 
obtained at a = 0.3 for n = l/v2, i.e. when the out- 
put temporal modes are produced with equal probabili- 
ties. Similar to the previous works [13, the Wigner 
function (22) can be experimentally reconstructed from 
the data of balanced homodyne detection, when the ui\ 
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FIG. 3: (color online) Plot of the combination B(J) (Eq.(32)) 
for three values of n: 0.4 (dashed); l/v2 (solid) and 0.9 
(dotted) corresponding to three values of £1 in Fig.2. (from 
top to bottom). 
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signals at the detectors are measured at two different 
times matched to the time separation between two out- 
put pulses obtained in Fig. 2. However, a rigorous experi- 
mental demonstration of maximal violation is attainable 
in a loophole-free Bell test, when both locality and de- 
tection loopholes are closed in a single experiment [HI]. 
The locality-loophole can be easily avoided, if the homo- 
dyne detection of two co-propagating temporal modes is 
performed by one-photon detectors as fast as the two si- 
multaneous measurements on each of the modes are sepa- 
rated by a spacelike interval. This may easily be achieved 
within our model, because the time separation between 
the two temporal modes amounts to several nanoseconds. 
To eliminate second, detection-efficiency loophole, we fol- 
low the works [13], although in these papers SPDC in 
nonlinear crystal has been used for generation of time- 
bin qubit. To satisfy the Banaszek-Wodkiewicz crite- 
rion one has to prepare pure single-photon state or the 
"vacuum-cleaned" Wigner function. To this end the dif- 
ferent methods have been applied in [27| leading to the 
same result and showing a strict violation of Bell's in- 
equality in accordance with theoretical expectations. 



IV. CONCLUSIONS 

In conclusion, we have demonstrated the possibility for 
dynamic preparation of a single photon in distinct tem- 
poral modes, employing strong parametric interaction 
between two slow single-photon pulses and their group 
delay. Disregarding the photon losses, we have found 
the solution of propagation equations for the quantum 
field operators depending on the propagation distance in 
terms of the Bessel function. Although the losses result in 
a decay of the field amplitudes, they do not prevent the 
temporal splitting of quantum pulses. Moreover, since 
the two well separated u>% pulses undergo the same losses, 
the entanglement between them is almost insensitive to 
losses and easier to purify, so that the proposed scheme 
can be regarded as a robust source of narrow-bandwidth 
single-photon qubits. We have shown the ability of our 
scheme to achieve an arbitrary entanglement between the 
toi temporal modes by adjusting the driving field inten- 
sity, while the separation between the time bins can be 
controlled by using the different atomic-level configura- 
tions to obtain the different group velocities of quantum 
fields. Subsequent papers will discuss a possibility for 
transferring and distributing quantum information be- 
tween optical modes of different frequencies ui\ and lj 2 in 
a loss and decoherence-free fashion, as well as the more 
complicated case of two input single-photon pulses will be 
analyzed and the results of detail numerical simulations 
will be presented. 

This work was supported by the ISTC Grant No.A- 
1095 and INTAS Project Ref.Nr 06-1000017-9234. 

APPENDIX A: COMMUTATION RELA- 
TIONS FOR FIELD OPERATORS 



In this Appendix we derive the commutation rela- 
tions Eq.(l) for traveling-wave electric field operators. 
In free space z < 0, the field operators are given by 

oo 

£i(z,t) = ai^ n e~ lUrl ^ t ~ z ^ c \ where di, n is the annihi- 

n=0 

lation operator for the i-th field discrete mode with a 
frequency uji +uj„. These operators satisfy the standard 
boson commutation relations 

[ai tn ,af n ,] = 5ijS nn > (33) 

In the continuum limit a^ n — > (Aw) 1 / 2 ai(w) and VJ — > 
(1/Aw) J du> with Alo — 2ttc/L, the commutation rela- 
tions at z=0 arc found to be 



[fi(o,t),ft(o,0] = -M(*-*0 



(Al) 



Inside the medium, the commutation relations satisfy the 
equation 

^[f 1 (z,i),£ 1 + (z,t')] = -^(| + | 7 )[fi(z,i),£+(^0] 



- ip([S 2 (z, t), £+{z, t')] - [&(z, t), £+(z, t')]) (A2) 

With the use of Eqs.(14) we have 

[£ 2 (z,t),£+(z,t / )} = [£ 1 (z,t),£+(z,t / )} = 

Recalling also that [£i(z 7 t),£^(z, t')] is a function of time 
difference t — t' and, hence, 

(| + | 7 )[£ 1 (z,t),£+(z,O]-0 

we obtain that the commutation relations are spatially 
invariant and have the form of Eq.(2). 

APPENDIX B: DESTRUCTIVE INTERFER- 
ENCE IN PULSE TEMPORAL SPLITTING 

The goal of this appendix is to show that the two pro- 
cesses responsible for lo\ pulse splitting interfere destruc- 
tively. For that we calculate the u)\ single photon wave- 
function 

$i(z,t) = (0 | | lr) = h(t-z/v 1 ) + 



+ / dxf 1 (t-z/v 2 x) — (Bl) 

J V1V2 OZ 



where Eq.(16) has been used. Using instead of a; a new 
variable y = x/z and taking into account 



dJo(i>) T i ..dtp I y 



where ip is now written as ip = 2(3z^Jy{\ — y), we finally 
get 

$i(z,t) = f 1 (t-z/v 1 )- 
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l 

-z/3 [ dyh{t-z/v 2 - V -l^zy)J 1 {^)>r^ (B2) 
J vxv 2 V 1 - V 



It is easy to see that ip varies within the limits < ip < (3z 
where, for the chosen parameters, (3z ^ pL ~ 3. Since 
in this interval Ji(ip) > 0, the integrand in the second 
term in Eq.(A2) is positive and, thus, the parametric 
regeneration of u)\ photon interferes destructively with 
the first term fi(t — z/v\), which describes the u)\ pulse 
slowing under the EIT in the ladder system. 

It is convenient to express $i(z, t) as a sum of two 
wavefunctions corresponding to the fast and slow im- 
pulses 

3>i(z,t) = $[(z,t) + $f(z,t) (B3) 



where $f' S (z, t) are real functions and $f (z, t) is always 
negative. By numerical calculations one can show that, 
for a given time, 

fltoiMWO (B4) 

and at a fixed z 

J $?[(z,t)$¥{z,t)dt = (B5) 

indicating that the two temporal modes are spatially and 
temporally well separated systems sharing one photon. 
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